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ABSTRACT 

Accretion onto a supermassive black hole of a rotating inflow is a particularly difficult problem to 
study because of the wide range of length scales involved. There have been broadly utilized analytic 
and numerical treatments of the global properties of accretion flows, but detailed numerical simula- 
tions are required to address certain critical aspects. We use the ZEUS code to run hydrodynamical 
simulations of rotating, axisymmctric accretion flows with Brcmsstrahlung cooling, considering solu- 
tions for which the centrifugal balance radius significantly exceeds the Schwarzschild radius, with and 
without viscous angular momentum transport. Infalling gas is followed from well beyond the Bondi 
radius down to the vicinity of the black hole. We produce a continuum of solutions with respect 
to the single parameter Ms/A/Edd, and there is a sharp transition between two general classes of 
solutions at an Eddington ratio of M-q/ M^d ~ few x 10~ 2 . Our high inflow solutions are very similar 
to the standard Shakura & Sunyaev (1973) results. But our low inflow results are to zeroth order 
the stationary Papaloizou and Pringle (1984) solution, which has no accretion. To next order in the 
small, assumed viscosity they show circulation, with disk and conical wind outflows almost balancing 
inflow. These solutions are characterized by hot, vertically extended disks, and net accretion proceeds 
at an extremely low rate, only of order a times the inflow rate. Our simulations have converged with 
respect to spatial resolution and temporal duration, and they do not depend strongly on our choice 
of boundary conditions. 

Subject headings: accretion: accretion disks - black hole physics - quasars: general - X-rays: general 



1. INTRODUCTION 

Supermassive Black Hole accretion is affected by the 
properties of the accretion flow on a wide range of length 
scales. Gravity and gas dynamics play an important role 
from the gravitational radius out to much larger dis- 
tances where the gravitational potential energy of the 
gas is comparable to its kinetic energy. The properties 
of the infalling gas at even larger radii, where gravity is 
weak can also in principle affect the global properties of 
the accretion flow, an important point that is sometimes 
overlooked in studies of accretion physics. 

In the unmagnetized spherically symmetric case with 
no angular momentum and adiabatic equation of state, 
the global ac cretion flow is given by the Bondi solution 
(jBondi 111953 ). The Bondi radius R B = GM B h/c^ j00 sets 
the scale at which gravity becomes important, and the 
mass accretion rate is given by Mb = A47rrgp 00 c Si00 . 
Mbh here is the black hole mass, is the mass density 
of gas at infinity, c Si00 is the gas sound speed at infin- 
ity, a nd A = [2/(5 - 3 7 )](5-3 7 )/[2( 7 -i)] / 4 flFrank et al.1 
120021) . where 7 is the adiabatic index. The introduc- 
tion of angular momentum fundamentally modifies this 
picture, however, and for sufficiently high angular mo- 
mentum the gas can no longer accrete onto the central 
black hole. For gas with specific angular momentum j, 
the centrifugal radius R c = j 2 /GMbh defines the ra- 
dius interior to which the gas has difficulty penetrating. 

jgli@astro.princcton.edu 



Accretion in this case requires some form of angular mo- 
mentum transport. Turbulent stresses associated with 
the nonlinear development of the MagnetoRotational In- 
stability (MRI) are tho ught to be capable of transport- 
ing angular momentum (IVelikhov I ll959t IChandrasekharl 
119601 : iBalbus fc Hawlevlll991|) . The angular momentum 
transport properties of the MRI are often parameterized 
with an effective shear viscosity v = ac s H p , where nu- 
merical simulations suggest a is typically of order 0.01 
dStone et al. II 19961: iHawlev et alJl2Cmt iMcKinnev et al.l 
l2012f ). Observations actually sugge st a vertica lly and 
azimuthally averaged a ~ 0.1 - 0.4 (|King ||2007| ). a dis- 
crepancy that to date has not yet been fully addressed. 

There are a number of general classes of analytic so- 
lutions for rotating flows with viscous angular momen- 
tum transport. The standard solution in the limit of a 
geometrically thin and optically thick viscous disk was 
given by Shakura & Sunyaev (1973; see also Lynden- 
Bell & Pringle 1974, Novikov & Thorne 1973), com- 
monly referred to as an a-disk. Thin disks typically 
neglect global heat transport and only consider effects 
that are first order in the disk thickness H/R, where H 
is the disk thickness at cylindrical radius R. The next 
order of approximation is slim disks, which account for 
terms in the equations of motion that are second order 
in the disk thickness H/R (Abramowicz et al. 1988; 
see also Abramowicz et al. 1986, Kato et al. 1988, 
Chen & Taam 1993, Narayan & Popham 1993, Katz 
1977, Begelman 1978, Begelman & Meier 1982, Eggum 
et al. 1988). Slim disks allow for non- negligible ra- 
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dial velocity, horizontal pressure gradients, and advective 
heat transport, and they have been used to model geo- 
metrically thick, hot disks accreting at super-Eddington 
rates. Advection Dominate d Accretion Flows (ADAFs , 



Naravan fc Y] 11994 Il995a|bl: lAbramowicz et all 119951: 
Ichimaru I Il977t iRees et al.l Il982t iQuataert fc Naravanl 



1999al INaravan fc Fabian" 120 111) describe radiatively in- 
efficient sub-Eddington flows. The accreting low den- 
sity gas is hot and geometrically extended in the po- 
lar direction, and the thermal energy is advected onto 
the central object. The closel y related Adiabatic In- 
flow Outflow Solution s (ADIOS, iBlandford fe Begelmanl 
ll99ft lBcgcrman 20TH) can accrete at much lower rates by 
allowing for outflows th at can carry off en ergy and angu- 
lar momentum (see also lJiao fc Wul[201lD . At present it 
is unclear which, if any, of these solutions real supcrmas- 
sive black hole accretion flows will pick. Certainly for hot 
disks with H ~ R, the standard approximation s in thin 
disk accretion theory break down (jPri nglc 198l|), and we 
should consider terms of all order in Hj R. Further study 
of the detailed properties of geometrically thick disks is 
required. 

Numerical simulations have provided further insight 
on accretion physics , and Convection Dominated Accre- 
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, which operate in the weak viscosity regime for 
flows with no cooling, have been proposed to ex- 
plain one c l ass of simulations. Recent work by 
lYuan et al.l (|2012bO tends to favor ADIOS-type mod- 
els over CDAFs, however. Other authors have stud- 
ied the geometry of two dimensional magnetohydrody- 
namical flows accreting f r om to ru s-like initial re s ervoir s 
of gas dStone fc Pringlel d200l: lHawlev et al.l (l200l[) : 
Machida et al.l (120011): llgumenshchev fc Naravanl (120021) 



Hawlev 



;t al.l llzUUll);llg umcnsn chev &i IN aravan ( zuua) : 
Balbusl (120021 ) : see also INaravan et al. I (|2012D 



for 3D simulations with general relativity). Many previ- 
ous works neglect the radiative properties of the flow, 
however, which may play an impor tant role in de- 
termining the geom etry of the flow (jDibi et al.1 120111 : 
iFragile et al. 1 120121) . For accretion flows with radia- 
tive cooling, the properties of the infalling gas at in- 
finity determine the strength of the radiative cooling 
and can affect the global properties of the flow (Yuan 
et al. 2000). Several authors m imic cooling by varying 
the adiabatic index of the ga s ( Moscib rodzka fc Progal 
120081 1200ft I Janiuket al ll2009D or include radiative tr ans- 
fer (lOhsuga et al. 11200a lOhsuga fc Mineshige 112011b see 
also lZanotti et al.ll2011 , and references therein), but ac- 
cretion in these works is aided by at least some inflowing 
gas that has low angular momentum or starting with 
gas that is bound to the black hole. Another important 
consequence of radiation from gas near the black hole is 
that it can drastically tra nsform the outer parts of th e 
accretion flow via heating (jPark k Ostriker Il200l [2007h . 

It is well known that observed black holes can accrete 
at rates that are only a small fraction of the Bondi ac- 
cretion rate. For example, Sgr A* appears to be ac- 
creting at a rate of 10~ 3 - 10" 2 Mb (|Yuan et al. 1 120031: 
IQuataert"! 12004) . A number of authors have studied 
the transition between the cooling flow at large radii 
and the accretion flow near and below the Bondi ra- 



dius ([QuataerF2!SiE^i3[illl3 IQuataert fc Naravan I 
120001 ). and the implications for accretion rates. Vari- 
ous energy transport mechanisms have been suggested 
that may be a ble to reduce accretion rates well below 
the Bondi rate dCiotti fc Ostriker [IpOOll: IQuataert 112004 
Sharma et al.l 120081 : iShcherbakov 1 120101 : llgumenshchev I 



20061 ). Other authors have invoked Radiatively Inef- 



ficient Accretion Flows to expl ain the low luminosity 
of ob served accreting systems ([Naravan fc McClintockl 
1200a and references therein). We consider here the al- 
ternate view that it is simply the angular momentum of 
the gas that slows accretion. There has been some ana- 
lytic and semi-analytic work explori ng this idea, includ- 
ing angular momentum tran sport ([Xu fc Chen I 119971 ; 
IBlandford fc Begelmanl I1999D . and both angular mo- 
mentum tra nsport and cooling (e.g.. iChakrabarti l l l996l : 
|Park1l27Mlino gamov fc Sunvaev N2010D. as well as nu- 
merical work without coo lin g (e.g.. | Stone et al. I IT999I: 
iProga fc Begelmanl l2003a|bl: lYuan et al. I l2012al H) and 
with cooling ([Yuan fc Bufl2010l)~ 

This paper will treat the accretion problem in the do- 
main where the angular momentum is not sufficient to 
prevent inflow at the Bondi rate but ample to inhibit 
accretion without the operation of some viscous process. 
That is, we will work exclusively in the domain 



R s < R c < R B , 



(1) 



where R s = 2GMbh/c 2 is the Schwarzschild radius. 
Since the ratio of the bounding radii, Rb/R s ~ 10 8 /T4, 
this choice is well-defined. Moreover, many ellipti- 
cal galaxies, including M87, have stellar populations 
with small average spin that is insufficient to cre- 
ate Keplerian disks near t he Bondi radius (see e.g., 
Ilnogamov fc Sunvaev 1 120 lot ). We run hydrodynamical 
simulations to follow infalling, radiating gas that starts 
from well beyond the Bondi radius with Be > c 2 00/(7 — 
1) (see Equation IT0|) down to the vicinity of the black 
hole. We include shear viscosity to capture qualitatively 
the angular momentum transport properties of the MRI, 
and so accretion can proceed at varying rates, depend- 
ing on the flow parameters. Since the gravitational and 
centrifugal forces do not have the same scaling with ra- 
dial distance in our setup, we do not expe ct our solutions 
to be self-similar as in sta ndard ADAF (INaravan fc Y 



1991 and CDAF models ([Quataert fc Gruzinovl 1200 



Naravan et al.ll20"o"oT ) . The thermal Bremsstrahlung cool- 



ing can also play an essential role in the energetics of 
the flow and determine if gas is bound to the black hole 
or not. We provide accretion flow solutions in both the 
very sub-Eddington, nearly adiabatic domain, and also 
solutions approaching the super-Eddington domain. Sec- 
tion [2] provides the mathematical framework framework 
for our problem and describes our numerical method. 
Section [3] gives our main results, and the broader impli- 
cations of our results are discussed in sections [4] and [5j 
The Appendix provides details on the convergence prop- 
erties of our solutions. We will designate the low inflow 
solutions that we find as "RRIOS" - Radiating Rotating 
Inflow-Outflow Solutions. 

2. METHOD 
2.1. Equations 
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We use ZEUS-2D v2.0 (Stone & Norman 1992) to solve 
the equations of hydrodynamics, 



d_ 

dt 



v-V)p + pV-v = 0, 



at 



p 



v • V ] v = - Vp - pV$, 



at 



v ■ V 



-pV • v, 



(2) 
(3) 
(4) 



where p is mass density, v is flow velocity, e is energy den- 
sity, gas pressure p = (7 — l)e, where 7 is the adiabatic 
index, and the gravitational potential $ = —GMbh/t. 
We run global hydrodynamical simulations of axisym- 
metric rotating accretion flows in spherical polar coordi- 
nates. The simulations include thermal Brcmsstrahlung 
cooling, shear viscosity to capture drag due to turbulent 
stresses, and modified boundary conditions to allow for 
disk and wind outflows. The source terms for momentum 
and energy in the ZEUS hydrocode are modified by the 
addition of the terms within square brackets below: 



dv 

' y~dt 

and 



- Vp - pV$ - V ■ Q + [V • a'] , (5) 



— pV • v — Q • Vi>- 

\<j' ■ Vl) - e B rcm] , 



(6) 



in which cr' is the viscous shear tensor and q is the heat 
flux due to thermal conduction of electrons. ZEUS' stan- 
dard tensor artificial viscosity Q is applied with shocks 
spread over m 2 zones, and we use a Foward Time Cen- 
tered Space differencing scheme for our modified source 
terms. We neglect the non-azimuthal components of the 
viscous shear tensor cr', as it is believed that the poloidal 
shear is subdominant to the az imuthal shear in th e non- 
linear regime of the MRI (see iStone et al. 1119991 ). The 
azimuthal components of the viscous shear tensor are 



/ 9 fvj, 

a rcj) = vpr— — 
or V r 



vpsm.6 d / 



86 Vsin# 



(7) 



(8) 



where v is the kinematic viscosity. The Bremsstrahhmg 
cooling term (Svensson 1982, Ball et al. 2001) 

e Br em = a f r 2 e m e c 3 n 2 (32/3) (2/tt) 1/2 (-^) 1/2 , (9) 

m e c^ 

where 07 is the fine structure constant and r e is the clas- 
sical electron radius. We use this abbreviated form for 
the Bremsstrahlung cooling to capture the qualitative 
effects of thermal emission; in fact in high temperature 
regions the electron component approaches relativistic 
speeds and Bremsstrahlung radiation increases substan- 
tially above equation ©. The cooling is limited by a 
floor in the temperature, equal to the gas temperature 
at the outer boundary, in order to maintain stability. The 
advective transport terms in ZEUS are unmodified. 



2.2. Setup & Flow Properties 

In our simulations gas flows in through the outer 
boundary, located at -Rout = IORb, at the Bondi rate, 
Mb • Wc modify the outer boundary condition locally to 
an outflow boundary condition in whichever grid zones at 
the outer boundary that the flow is moving radially out- 
ward, however. We note that several other authors have 
studied accretion flows with outer boundary near th e 
Bondi radius fe.g. JPen et al. 1 120031: iPang et al. 112011ft . 
Our inner boundary, located at R ln = 10 _3 i?s, is a 
standard no torque outflow boundary condition. i?j n is 
located at 91R S , and we have checked the convergence 
of our results with respect to the location of the in- 
ner boundary (see Appendix). The flow is initialized 
throughout with Bondi profile for density, energy den- 
sity, and radial velocity. Our standard simulations are 
run with = 2 x 10 7 K. This is close to the tempera- 
ture of gas heated radiatively by a typical quasar output 
spectrum ( Sazonov et al. 1 120050 . The simulations are 
run with adiabatic index 7 = 1.5, except where other- 
wise noted, and we verify explicitly that this assump- 
tion does not affect our main results (cf Figure [7]). The 
real gas will have 7 = 5/3 above R ~ 300 — 500i? s . At 
lower radii, in the non-collisional case, electrons become 
relativistic while protons are nonrclativistic, leading to 
lower effective 7. Our choice for 7 is made in order 
to allow us to well separate the sonic and centrifugal 
radii for supersonic infall. The density at infinity, p^ , is 
a parameter that determines the dimensionless quantity 
Me/MEdd; where MEdd = AirGM^Rm p /eo-TC is the Ed- 
dington luminosity and we assume an efficiency e = 0.1. 
This dimensionless mass accretion rate parameter scales 
as Me/MEdd oc MbhPod/cI ^. The dimensionless cool- 
ing fraction e B rcmAi char /nfcT cx M B hAx>/Cs )00 , where we 
have taken characteristic time Ai c h ar = M-qh/Cs qo- The 
black hole mass and gas density at infinity appear to- 
gether in the above two quantities, and the free parame- 
ter in the accretion flow is actually the product MbhPoo 
(see Chang & Ostrikcr 1985). This parameter effectively 
determines the strength of the Bremsstrahlung cooling, 
with Ms/A/Edd <§C 1 corresponding to no cooling and 
Mb /-MEdd S> 1 corresponding to strong cooling. For 
convenience we pick Mbh = 2 x 10 8 Af Q and allow p^ to 
vary freely. At p^ = 2.3 x 10~ 21 .g cm -3 , the Eddington 

ratio Afe/MEdd = 0.1. 

Gas is initialized with constant specific angular mo- 
mentum throughout, except near the rotation axis, where 
it tapers to zero. The inflowing gas at the outer bound- 
ary has the same angular momentum, without tapering, 
ensuring that after initial transients, the gas has uniform 
angular momentum throughout. For adiabatic inviscid 
flows, the Bernoulli constant 



Be = \ (vl 



ill 
2w 2 



GM 



7-1 



(10) 



where w is cylindrical radius, is a conserved quantity 
(see e.g., Narayan & Yi 1994). It immediately follows 
that there is a centrifugal radius interior to which the 
gas cannot penetrate if gas starts from rest at infinity 
with Be = 7c 2 ,00/(7 — 1). The specific angular momen- 
tum is another free parameter in our simulations, but 
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we pick joo = v / 0.02i?BC s ,oo so that the centrifugal ra- 
dius R c — 0.02i?B is well resolved and lies at a radius 20 
times that of the inner boundary. Hence we are comput- 
ing within the domain given by equation ([1]). The limit 
of weak angular momentum with R c -C R ln approaches 
the Bondi solution. For our choice of specific angular mo- 
mentum, the introduction of shear viscosity and angular 
momentum transport allows for gas to migrate from the 
centrifugal radius down to the inner boundary and to 
accrete. The accreting gas carries only a small fraction 
(< 20%) of the initial angular momentum through the 
inner boundary. In our simulations we set the kinematic 
viscosity v to a constant nonzero value v = 10 _3 c SiOO i?B- 
The effective a coefficient, a = v/c s H p , where c s is the 
local sound speed and H p the local pressure scale height 
(see Section l3.3|) . is then spatially variable, but is typi- 
cally of order ~ 0.0 1 near the centrifug al barrier where 
the density peaks. IStone et al. I (jl999f ) found that the 
properties of their non-radiative hydrodynamic flow do 
not depend s trongly on the as s umed form for kinematic 
viscosity, and IStone fc Pringld ([200l found weakly mag- 
netized flows to be qualitatively similar to the viscous 
hydrodynamic ones. Incidentally, the limit of large spe- 
cific angular momentum with R c > i?e should approach 
that of a giant viscous disk with large radial extent, but 
we have not explored this regime in any detail. The 
Bremsstrahlung cooling also modifies the flow profile, 
and for cold disks the inflow can be highly supersonic. 
In these cold solutions the sonic radius is typically inter- 
mediate between but still well separated from the cen- 
trifugal and Bondi radii. 

2.3. Numerics 

Our standard simulations are run with 21 equally 
spaced grid zones in polar angle 9 varying from 9 = 
to 9 = 7r/2, enforcing symmetry across the equatorial 
plane. In the radial direction we use a non-uniform log- 
arithmic grid with 32 zones per decade ranging from 
i?in = 10~ 3 i?B to i? ut = 10i?s, with 128 total grid zones. 
As a code check, we ran the spherically symmetric test 
problem with no cooling. Our code is able to maintain 
the Bondi profile for mass density and radial velocity to 
within 5% for various values of the adiabatic index 7. 
The internal energy density is less well conserved in the 
very inner portions of the flow (« 10 innermost cells) 
due to artificial viscous heating, but is still within a fac- 
tor of ps 2 of the Bondi solution. The simulations are 
run until all time-averaged gross quantifiers of the flow 
structure approach quasi-steady values. All diagnostics 
of the flow are averaged over many tens of Bondi times 
tB = Rb/cs.oo, typically between times t = 40 £b and 
t = 90 tB (see Appendix). Flow patterns are illustrated 
at time t = 90 t-g, except where otherwise noted. The 
code is mass conservative, so mass flowing inwards from 
large radii must accumulate on the grid or outflow either 
through the inner boundary or back out to large radii via 
the combination of polar and equatorial outflow. 

We also run a number of tests to check the proper im- 
plementation of the additional physics we have added to 
ZEUS. We verified the proper thermal cooling of a sta- 
tionary constant density g well as the diffusion of 
thermal energy when a temperature gradient exists. Fi- 
nally, we have reproduced qualitati vely the viscou s evo- 
lution of an initial ring of matter (jPringle I HMD) . The 



ring spreads in radius as matter moves inwards, but the 
bulk of the angular momentum is transferred outwards. 
Kinetic energy lost to viscous dissipation is added to the 
internal energy of the gas, with the viscous substep of 
the code conserving energy to within 10%. 

3. RESULTS 
3.1. Adiabatic Inviscid Flows 

For adiabatic flows with sufficiently high angular mo- 
mentum that the centrifugal barrier is located far from 
the black hole, equation (fpQ|) immediately tells us that 
there can be no accretion onto the black hole. Any inflow- 
ing high angular momentum material must either accu- 
mulate beyond the centrifugal radius or be deflected back 
out wards and ejected in outflows (iProga fc Begelnia " 
see | 20 3a|. and references therei n; see also lHawlev et 
Il984allbl : lHawlev fc Smarrl [19861) . 

In the adiabatic regime with constant angular momen- 
tum j throughout, there is a stationary (v r = vg = 
0) s ettling solution satisfying the momentum equation 
(see Papaloizou k, Pringle 1 119841 : iFishbone &: Moncrief I 
11976 ). Starting from the momentum equation, 



1 f 
Vp - V$ + ^tSj = 0, 

p w 6 



the adiabatic relation p cx p 1 gives 



3 



= 0. 



k 7- 1 ' ' 2ro 2 ; 
This yields an equation for the sound speed, 



. GM 
1+ 7-1-2 



7-1 3 

2 C* TJ7 3 



(11) 



(12) 



(13) 



valid where the right hand side is positive, as well as the 
density profile, 

f = [l + (7-l)^;-V^ 1/(7 " ) - < 14 > 
Poo cj^r 2 cj^-cu 2 

The density is zero near the polar axis in the region 
where the right hand side of equation (|T3| is nega- 
tive. This solution is actually d ynamically unstable to 
non-a xisymmetric perturbations (jPapaloizou fc Pringle I 
fl98l . which we do not allow for in our axisym- 
metric code. Neglecting these instabilities is likely 
a reasonable assumption, however, as we have angu- 
lar momentum gradients th at can stabilize the flow 
(jPapaloizou fc Pringle 1I1987D . 

Our multidimensional adiabatic simulations closely re- 
semble this stationary settling solution, with sound speed 
profile matching equation (|13[) to within a few percent 
exterior to the centrifugal barrier. Figure la shows the 
structure of the flow within the Bondi radius at Edding- 
ton ratio Ms/MEdd = 10~ 4 . The red contours are log- 
arithmic contours of density, given in cgs units, and we 
overlay black contours of logarithmic density for the adi- 
abatic settling solution. The solid black line shows the 
zero density surface in the stationary solution. Blue ar- 
rows indicate vectors of mass flux in the R— z plane, with 
length proportional to the magnitude of the mass flux 
vector. Subsequent diagrams of flow patterns showing 
logarithmic density contours and mass flux vectors are 
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illustrated in a similar fashion. The pressure associated 
with the accumulation of matter in an extended torus im- 
pedes the inflow of matter, as in the stationary solution. 
Indeed, the mass inflow rate at the Bondi radius is re- 
duced by more than one order of magnitude to ~ 0.06Mb 
at t = 45£b- There is evidence of time-dependent circu- 
lation, however, partially due to our choice of initial con- 
ditions, and the flow is not purely stationary with zero 
inflow rate at the Bondi radius. Matter flows inwards 
at intermediate latitudes and turns around at radii ex- 
terior to the centrifugal barrier, flowing back outwards 
in the polar and equatorial regions. Figure lb shows the 
flow structure interior to 0.1i?B, with centrifugal radius 
located at R c = 0.02i?B- There is again evidence of time- 
dependent circulation, but the largest differences in the 
flow relative to the stationary solution occur near the 
polar axis where the stationary solution has zero mass 
density. There are very small inflows in the polar region 
due to angular momentum transport associated with nu- 
merical viscosity, but the accretion rate is extremely low, 

-^accretion < 10 _3 Mb, i-e. any accretion is due to numer- 
ical error. In the absence of physical angular momentum 
transport the gas in the centrifugal barrier is stationary 
in the r and 9 directions, cannot accrete, and its pressure 
prevents new gas from falling down to small radii. Our 
results indicate that in the limit of very long and accurate 
integrations the solution would approach the stationary 
solution given in equations (|13|) and (fT4)) . since there is 
no engine to drive circulation for adiabatic flows. 

3.2. Hot and Cold Disks 

We find two general states for rotating axisymmctric 
accretion flows with radiative cooling and viscous angular 
momentum transport: hot thick disks with low accretion 
having polar and disk outflows, and cold thin disks with 
high accretion and weak polar and disk outflows. 

Figure 2a shows the flow structure interior to the Bondi 
radius for the hot disk with Mr/Mrm = 10~ 4 . Red 
contours are logarithmic contours of density, blue vec- 
tors represent mass flux in the R — z plane, and the solid 
black line indicates the zero density surface in the adi- 
abatic settling solution. Gas flows inwards over a wide 
range of intermediate polar angles, and there are equa- 
torial and polar outflows. The mass ejected in the polar 
outflow, ~ 0.02Mb, is similar in magnitude to the adia- 
batic case, and is again caused by inflowing matter turn- 
ing around exterior to the centrifugal torus. The detailed 
properties of this circulation may depend on the length 
and accuracy of our numerical integrations. The equa- 
torial outflow is driven primarily by azimuthal stresses 
associated with the shear viscosity, and it has magni- 
tude ~ 0.25A/b at the Bondi radius. The mass inflow 
rate through the Bondi radius increases over that in the 
inviscid case and essentially balances this equatorial out- 
flow. Figure 2b shows red contours of specific angular 

momentum in units of joo = V0.02i?BC s .oo, with blue 
vectors indicating the flux of angular momentum. The 
purple contour shows where j = joo] equatorward the 
gas has higher specific angular momentum, and poleward 
the gas has lower specific angular momentum. The equa- 
torial outflow carries of f the angular momentum of the 
gas that was accreted (iKolvkhalov fc Svunvaev 1 119791 : 
llnogamov fe Sunvaev 112010ft and some of the angular mo- 



mentum lost by the polar outflow. 

Figure 3a shows logarithmic density contours and ve- 
locity vectors interior to 0.1i?B- The black line indicates 
the zero density surface in the stationary solution. In- 
flowing matter dumps angular momentum into the equa- 
torial regions, and an outflow is driven beginning just 
beyond the centrifugal radius, located at R c = 0.02i?B- 
The inflowing gas is able to penetrate below the centrifu- 
gal radius only due to the angular momentum transport, 
i.e. the assumed weak MRI induced viscosity. In Fig- 
ure 3b, which shows logarithmic density contours and 
velocity vectors interior to half the centrifugal radius, we 
see that the gas accretes over a range of latitudes from a 
time-dependent, vertically extended, sub-Keplerian disk. 
The accretion rate is very low, however, with magnitude 
~ 4 x 10~ 3 A/b of order a times the mass inflow rate. 
There is some accretion due to the angular momentum 
transport, but the bulk of the gas is too energetic, with 
Be > c 2 s 00/(7 — 1), and has too much angular momen- 
tum for much of it to a ccrete. This conclus ion is in close 
agreement with that of IStone et al. I (|1999f) . 

In the cold disk solutions, the non-adiabatic gas radi- 
ates energy and most of the gas cannot travel back out to 
infinity on energetic grounds since Be < c 2 s 00/(7— 1). In 
the absence of angular momentum transport the gas con- 
tinuously accumulates near the centrifugal barrier. Gas 
pressure is insufficient to drive matter inward onto the 
central black hole and there is no steady solution. The 
introduction of angular momentum transport allows for 
significant accretion from the centrifugal torus, however. 
Figure 4a shows logarithmic density contours and mass 
flux vectors for the cold disk solution at Eddington ratio 
Ms/A/Edd =0.1, interior to the Bondi radius. The solid 
black line shows the zero density surface in the adiabatic 
settling solution given in equation (| 14|) . Matter accu- 
mulates in an extended torus whose pressure slows mass 
inflow, similar in spirit to the adiabatic settling solution. 
Inflow proceeds at low polar angles 8 < ir/4, however, 
with total inflow rate ~ 0.45Mb at the Bondi radius. 
Equatorial gas exterior to the centrifugal barrier is still 
driven outwards by the viscous engine, but the outward 
propagation of the equatorial outflow is slowed by the 
strong radiative cooling. In the absence of external pres- 
sure, the equatorial out flow would continue to propagat e 
outwards indefinitely (jKolvkhalov fc Svunvaev I 119801) . 
Our constant pressure outer boundary condition slows 
the outflow, however, and we find that the equatorial 
outflow is unable to propagate to the Bondi radius on 
time scales of order ~ 100 is- This solution is not tech- 
nically a steady-state solution, as angular momentum is 
continuously dumped into the radially extended torus, 
which acts as a large sink of angular momentum. If we 
integrated for much longer times, however, the equatorial 
outflow must continue to propagate outwards to conserve 
momentum in a steady-state. Figure 4b shows logarith- 
mic density contours and mass fluxes interior to 0.1i?B- 
The flow forms a geometrically thin disk, and angular 
momentum transport allows matter to continue to flow 
inwards from this disk. We illustrate the cold thin disk 
near the centrifugal barrier in greater detail in Figure 
5. Red contours are again logarithmic density contours 
and blue vectors represent mass flux. Accretion proceeds 
through the inner boundary from the cold, thin, sub- 
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Fig. 1. — (a) The structure of the inviscid flow with no cooling, interior to i?e at t = 45tB- Red contours arc logarithmic contours 
of density , an d blue arrows indicate mass flux in the R — z plane. The flow approaches a stationary settling solution (black contours, cf 
equation 1141 1 in which the inflow rate through the Bondi radius is reduced by an order of magnitude below the Bondi rate. The solid 
black line shows the zero density surface in the stationary solution, (b) Flow structure interior to 5R C = 0.1-Rb- The bulk of the gas in the 
disk cannot accrete and the accretion through the polar funnel is negligible and due to numerical error. 
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Fig. 2. — (a) The structure of the hot disk with angular momentum transport designed to mimic the MRI, shown interior to i?B tn the low 
accretion limit where Ma/MEdd = 10 -4 . Cooling is negligible. Red contours are logarithmic contours of density, blue arrows indicate mass 
flux in the R — z plane, and the black line demarcates the zero density surface in the adiabatic settling solution. Matter flows in through the 
Bondi radius and back out in disk and conical polar outflows, (b) Contours of specific angular momentum in units of joo = \/0.02_Rb c s,oo 
and blue vectors of angular momentum flux. The black line demarcates the zero density surface in the stationary solution, and the purple 
contour marks where the specific angular momentum j = jao. Angular momentum of inflowing material is transferred to the equatorial 
outflow. 
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Fig. 3. — Closeup of the structure of the hot disk with angular momentum transport at Eddington ratio M-q/ M-^dd = 10 -4 . This is the 
same model as Figure 2. Red contours are logarithmic contours of density, and blue arrows represent mass flux in the R — z plane, (a) Flow 
structure interior to 5R C = 0.1-Rb- The black line shows the zero density surface in the adiabatic settling solution. Angular momentum 
transport allows a small amount of matter to inflow below the centrifugal radius, and viscous stresses drive a strong equatorial outflow, 
(b) Flow structure interior to half the centrifugal radius or 0.01-Rb- Accretion proceeds at very low rates over a range of polar angles. 
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Keplerian disk, at a rate nearly equal to the mass inflow 
rate through the B ondi radius. This case res embles the 
standard thin disk (Shakur a fc Sunvaevlll973l ). 

3.3. Radiative Cooling 

Figure |6] shows the quasi-steady mass accretion rate 
through the inner boundary, the equatorial and polar 
wind mass outflow rates through the Bondi radius, and 
the mass inflow rate through the Bondi radius, all as 
functions of the parameter Mb / M^dd- Crosses on the ac- 
cretion rate curve indicate the Eddington ratios for which 
our simulations have converged with respect to tempo- 
ral duration. For Me/A/Edd < 0.01, the Bremsstrahlung 
cooling is weak and plays little dynamic role in the flow. 
The flow converges to a low accretion rate solution for 
Me/MEdd ^ 10~ 3 . The non-zero viscosity allows only 
a small fraction, ~ 0.4%, of the inflowing matter at the 
Bondi radius to accrete onto the central black hole. The 
mass inflow rate at the Bondi radius Mi n a ow is itself only 
~ 0.29Mb because the centrifugal torus and outflows sig- 
nificantly impede the inflow (cf Figures 2 and 3). The 
gas density on the equator is much larger than in the po- 
lar region, and the equatorial outflow dominates the out- 
ward mass and momentum flux at the Bondi radius. The 
equatorial outflow nearly matches the mass inflow rate 
at the Bondi radius, whereas the polar outflow is much 
less significant and an order of magnitude weaker than 
either. It is important to note, however, that viscous 
heating near the centrifugal radius leads to a buildup 
of thermal energy and can lead to intermittent episodes 
of convective overturning. The overturning episodes can 
last for ~ 10ts, with outward mass fluxes over a broad 
range of polar angles and a significant fraction of the 
characteristic mass inflow rates. The filled boxes in Fig- 
ure |5] represent estimates of the time-averaged polar mass 
outflows including these overturning episodes, which by 
definition are not included in our quasi-steady average 
mass fluxes (see Appendix for more detail). 

The strength of the radiative cooling increases with 

the parameter Me/MEdd, and in the range Ms/MEdd ~ 
0.01 — 0.1 there is an abrupt transition in the prop- 
erties of the flow to the strong cooling regime. For 
Me/MEdd ~ 0.1 the inflowing gas is impeded in the 
equatorial regions by gas pressure in an extended cen- 
trifugal torus (cf Figure 4) , and matter flows in through 
the Bondi radius at roughly half the Bondi rate. Nearly 
the entirety of this inflowing gas is accreted onto the 
centr al black hole and the solut ions approach the stan- 
dard iShakura fc Sunvaevl (|1973l ) picture. Most of the 
inflowing gas does not have enough energy to travel back 
out to infinity, and only a very small portion bounces 
off the centrifugal barrier and is ejected in a polar wind. 
The persistence of the polar wind is uncertain, and much 
higher resolution studies with additional physics such 
as magnetic fields and radiation are necessary to study 
the detailed properties of polar outflows. We empha- 
size again that the high Eddington ratio solutions are 
only quasi-steady (see Section l3.2j) . and the equatorial 
outflows would be stronger if we integrated for times 
> 1000i? B /c s ,oo. 

We have further verified that the low accretion rate 
solutions depend weakly on our choice of adiabatic index. 
Figure [7] shows the mass inflow rate from large radii, 



equatorial and polar outflow rates at i?B, and accretion 
rate as a function of 7. Mass fluxes depend weakly on 7 
in the range 1.4 — 1.65, and our ansatz 7 = 1.5 should 
not significantly affect our results. 

The high accretion rate flows from Figure [6] with 
Ms/MEdd > 0.1 are cold and form a geometrically thin 
disk, and the low accretion rate flows with Mb / Afedd < 
0.01 are hot and form vertically extended disks. Fig- 
ure 8a shows the emission weighted temperature, 



< T >= - 



,TdV 



J eBromdV 



(15) 



as a function of M-q/ M^dd- We show three characteristic 
temperature curves, from top to bottom averaged over 
the region interior to the centrifugal barrier, the total 
region within 0.1i?s, and in the centrifugal torus. The 
centrifugal torus is extremely dense and cools to a lower 
temperature (the floor that we set) than other parts of 
the flow at high Eddington ratios. The gas interior to the 
centrifugal barrier is hottest because the gas has pene- 
trated deeper into the gravitational potential well. For 
low Me/A^Edd < 0.01 the gas in the inner parts of the 
flow is hot and a sizable fraction of the centrifugal tem- 
perature 



T = - £ — — 100T 

hi R c 



(16) 



where fi = 0.5 is the mean molecular weight. Elec- 
trons become weakly relativistic at such temperatures, 
but we neglect in our computations the increase of the 
Bremsstrahlung energy losses for relativistic electrons. 
We measure the pressure scale height of the disk as the 
vertical distance at the centrifugal radius over which the 
pressure drops by a factor of e. The pressure scale height 
of the disk for the hot solutions with vertically extended 
disks is H p /R c ~ 0.75 (see Figure 8b), and they have 
a ~ 0.01 at the centrifugal barrier. The high accre- 
tion rate solutions with cold thin disks have pressure 
scale at the centrifugal radius H p /R c ~ 0.15, covering 
~ 5 cells in the highest resolution simulation (see Ap- 
pendix). For a cold thin disk in hydrostatic equilibrium 
and with angular momentum j = \/0.02i?BC SiOO we ex- 
pect Hp/ R c = Cs/vrf, = 0.07. Since we fix the kinematic 
viscosity v in our code, our cold disks have a ~ 0.4 rather 
than a ~ 0.01 at the centrifugal barrier. The code has 
difficulty running the cold solutions at lower viscosity, 
as there is significantly greater mass deposited near the 
centrifugal radius before accretion can proceed. The gas 
has cooled strongly and energetically cannot reach infin- 
ity, however, and we expect it to accrete independently 
of the strength of the viscosity, as long as it is non-zero. 
In any case, radiative transfer becomes important at Ed- 
dington ratios of ~ few x 0.01, physics that we do not 
include here. The precise location of the transition be- 
tween thin and thick disks, as well as the density, lumi- 
nosity, and optical depth of the thin disks, will depend 
on these choices. 

The conical polar outflows, driven by inflowing mat- 
ter bouncing off the centrifugal barrier, are in general 
weak, but they exist over a wide range of Eddington ra- 
tios. Figure 8c gives the half-opening angle of the polar 
outflow at the Bondi radius as a function of Me/A/Edd- 



10 




Fig. 4. — (a) The structure of the cold disk with angular momentum transport interior to _Rb, at Eddington ratio / M Edd = 0.1. Red 
contours are logarithmic contours of density, and blue arrows indicate mass flux in the R — z plane. Inflow is impeded by gas pressure at 
large polar angles but proceeds through the Bondi radius at polar angles 9 < tt/4. Angular momentum is dumped into the extended torus, 
(b) Flow structure interior to 5R C = 0.1-Rb- Matter continues to flow inwards from a geometrically thin disk. This is the case which most 
closely approximates the standard Shakura-Sunyaev disk. 
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Fig. 5. — A closeup view of the high inflow case with Me/A^Edd = 0.1. We show logarithmic density contours and mass flux vectors 
for the c old disk with angula r mom entum transport, interior to half the centrifugal radius. Matter accretes from the thin disk as in the 
standard Shakura & Sunyacv (1973) picture with angular momentum transport due to the MRI. 
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Fig. 6. — Mass inflow rate from large radii, equatorial outflow rate, polar outflow rate (boxes are averages that include episodes of 
convective overturning) , and accretion rate as a function of the Eddington parameter Mb / Mgdd • At high Mb / Mgdd mass inflows at ~ 1 /2 
the Bondi rate and nearly all of it is accreted. At low Ma/MEddi mass again inflows at ~ 1/2 Bondi, but nearly all of this matter is 
deflected outwards in disk outflows. There is an abrupt transition between these distinct solutions in the range M-g/M^dd ~ 0.01 — 0.1. 



The half-opening angle is typically 6*h a if-opcn ~ 5°. Wc 
show also the half-angles at which half of the mass and 
half of the momentum in the polar outflows are ejected, 
again computed at the Bondi radius. The zero density 
surface in the adiabatic settling solution lies at an an- 
gle of ~ 3° at the Bondi radius. Figure 8d shows the 
mass-flux weighted average wind velocity in the polar 
outflow, measured at the Bondi radius, as a function of 
ilie/A/Edd- The wind velocity can be a significant frac- 



tion of the centrifugal velocity 

v c = (2GM B n/Rc) 1/2 
with typical wind velocity < 



(17) 



n H >~ 0.15 v r . Thi s 
value is similar to that found in (|Yuan et al. Il2012bh . 
The outflow opening angle and velocity depend on the ra- 
dius at which we measure these quantities. The outflow 
propagates into a constant density and pressure medium 
in our setup, and at sufficiently large distances the out- 
flow will transfer its momentum and energy to the ex- 
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Fig. 7. — Mass inflow rate from large radii, equatorial and polar outflow rates at Rb, and accretion rate as a function of adiabatic index 
7. All solutions are computed for fixed Eddington ratio M^/M^dd = 10 — 4 , the hot disk solutions. Mass fluxes have relatively weak 
dependence on adiabatic index in the range illustrated, and we set 7 = 1.5 throughout the rest of this work. 



ternal medium. Beyond this region the outflow would 
become subsonic and probably circulate back to join the 
inflowing gas at intermediate latitudes. We quote values 
at the Bondi radius to give a quantitative estimate for 
the outflow properties at a definite radius. 

We now discuss in greater detail the radiative prop- 
erties of our accretion flows. If a spherically symmetric 
distribution of gas is in free-fall with velocity 



v s = (2GM B hA) 1/2 
and number density 



n(r) 



(r/Rs) 



-3/2 



(18) 



(19) 



where rj is the mean mass per electron, then the optical 
depth 



a T / n(r)dr cx M B /M 



Edd- 



(20) 



We compute the optical depth along three different sight- 
lines with fixed polar angle 9 = 0°,45°,90°, starting 
from the inner boundary at i?j n = 91i? s . We find that r 
scales linearly with Mb /M^dd (see Figure 8e) . The opti- 
cal depths approach 10 -2 just below M^/M^dd = 1CT 2 , 
where the flow switches between the hot, vertically ex- 
tended disk and the cold, thin disk, and there is a jump 
by a factor of ~ 100 in T90 at this transition. Further, 
making use of the temperature profile of an adiabatic gas 
in free-fall, 



T~T{R s )(r/R s ) 



-3( 7 -l)/2 



the total integrated Bremsstrahlung luminosity 



-^Brcm/ -^Edd OC 



1 



n ( r )2 T l/2 r 2 dr ( 22 ) 



^Edd 

cx (Af B /M E dd) 2 . (23) 

Figure 8f shows the Bremsstrahlung luminosity in our 
simulations integrated over the regions between R m = 
91i? s and 0.1i?B, and between R- ln = 91R S and the cen- 
trifugal torus. We obtain the scaling iBrem/^Edd oc 
(M-q / M-Edd) 2 for both integrated luminosities, and 

LBrcm — > 10" 7 iEdd at Mb / Mftdd = 10~ 2 . As we move 
the inner boundary at R[ n closer to the black hole, we 
expect the luminosity, temperature, and optical depth 
interior to the centrifugal radius to increase. Other work 
including general relativistic and magnetohydrodynamic 
physics and extending down to the horizon for non- 
spinning black holes suggests that essentially all matter 
flowing inwards at our inner boundary R; n = 91 R s will 
ultimately accrete, however (|Naravan et al. II2012D . 

4. DISCUSSION 

Observations suggest that the Active Galactic Nu- 
clei (AGN) in nature form a bimodal distribution, 
with broad-line AGN having accretion rates ~ 2 or- 
ders o f magnitude higher than narrow- line and lineless 
AGN (ITrump et al. 1 12001 [2ijTl IKollmeier etl!~ll200l 
iRussell et al.ll2012l ). The transition between these classes 
occurs at an Eddington ratio between 0.01 and 0.1, in 
remarkable agreement with our transition between hot 
and cold flows. Further, the luminosities of the cold 
and hot disk states suggests a natural feedback loop 
that can result in switching between the two solutions. 
From Figure [8] we see that the Bremsstrahlung emis- 
sion from a cold thin disk can approach 10 _2 LEdd for 
(21) Ms/M^dd ~ 0.1. Wc do not include the radiative fore- 
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Fig. 8. — 6 characteristic quantities in our solutions, all as a function of M^/M^dd- ( a ) Emission weighted temperature, averaged over 
the region interior to the centrifugal torus, the region interior to 0.1-Rbi and the centrifugal torus. The low accretion rate solutions are hot, 
with temperature a sizable fraction of the centrifugal temperature, whereas the high accretion rate solutions are cold, (b) Pressure scale 
height. The hot disks are vertically extended, with H v /R c ~ 0.75, whereas the cold disks are geometrically thin, with H p /R c <C 1. (c) 
Polar wind half-opening angle, measured at Rg, typically of order 5°. We show also the half-angles at which half the mass and momentum 
in the wind are ejected, (d) Mass-flux weighted polar wind velocity, measured at R^, typically a sizable fraction of the centrifugal velocity, 
(c) Optical depth along three sightlines with fixed 9 = 0°, 45°, 90°. r scales linearly with Mg / M^ddi but there is a jump in rgrj when the 
disk switches to a dense, cold thin disk. The optical depths approach 10~ 2 at Mb / -Mgdd = 10~ 2 . (f) Total Bremsstrahlung luminosity 
integrated over the region between the inner boundary and 0.1-Rbi and between the inner boundary and the centrifugal torus. Luminosity 
scales as the second power of M^/M^dd an d approaches 10~ 7 L^dd at A^B/A^Edd = 10~ 2 . 
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ing or thermal heating effects of this emission on in- 
falling gas in our code (park fc Ostriker 1 120011. 120071: 
I Yuan et al. 1 120091 ), but the Bremsstrahlung luminosity 
would significantly heat the inflowing matter. The disk 
could transition to the hot, vertically extended state with 
low accretion, effectively shutting off the Bremsstrahlung 
emission. In the absence of irradiating photons, the hot 
disk could then cool and begin accreting at high rates 
again, thus forming a cycle. This process may be an 
important detail in exp laining the high and low accre- 
tion states of quasars (jCiotti fc OstrikeTl 120011 ) . The 
physics of radiative feedback on supermassive black hole 
accr etion flows has been s t udied numerically in some de- 
tail (Sazonov et al. I l2005t lProgall2007t iNagamine et all 
201li iCiotti et al.ll2009l 120101: Kurosawa fcProgal 120091: 
Proga et al I" l200l . and a potentially significant effect 
here would be to increase the strength of the polar out- 
flows. 

Independently of any radiative feedback loops, our re- 
sults indicate that temporal variations in the external 
boundary conditions at infinity would lead to strongly 
time-dependent accretion flows. Observations are now 
approaching sufficient accuracy that for nearby super- 
massive black holes both the external boundary condi- 
tions surrounding the black hole and the reaction of the 
flow to these ex ternal boundary condit ions can be mea- 
sured. Indeed, iGillessen et al. I (I2012D recently discov- 
ered a gas cloud of only several Earth masses orbiting Sgr 
A* and with radius of closest approach ~ 3100i? s . The 
computations provided here are detailed enough to make 
concrete predictions concerning the geometrical, dynam- 
ical, and thermal properties of accretion disks and out- 
flows of nearby supermassive black holes from currently 
feasible observations. We do not provide detailed pre- 
dictions here but may further explore these issues in a 
separate paper addressing observational consequences. 

The lack of magn e tic fields is a limitatio n of our work 
dPessah et all 120081: iBeckwith et all 120091 iPenna et all 
2010t iTchekhovskov et all 1201 It iMcKinnev et alll2012D . 
Proga fe Begelmanl (|2003bl ) considered the non-radiative 
two-dimensional magnctohydrodynamic accretion prob- 
lem, and they find angular momentum transport and ac- 
cretion driven by the MRI. The two-dimensional struc- 
ture of the flow determines the accretion properties, with 
polar outflows quenching the accretion rate. Similar 
qualitative behavi or is present in the Adiabat ic Inflow 
Outflow Solutions (iBlandford fe Begelmanll999| ) , and we 
also see circulation and polar outflows in our hydrody- 
namical simulations with weak cooling. In principle the 
nonlinear properties of the MRI could modify the de- 
tailed structure of the flow in the weak cooling regime. 
As the cooling rate increases, however, the total energy 
of the fluid determines its fate and the detailed kinetic 
properties are less important. 

Another extension of our work relates to the en- 
forcement of reflection symmetry through the equatorial 
plane. The flow does not necessarily need to preserve 
this symmetry, and in principle there could be asymmet- 
ric, unstable modes that seed turbulence and/or convec- 
tion that modifies the properties of the accretion flow. 
We have explored this idea in some detail, and we find 



that the main difference from the presented results to 
be that the strong equatorial outflow of high angular 
momentum gas n eed no longer be ce ntere d exactly on 
the equ ator, as in iStone et al. I (|1999f ) and I Yuan et al. I 
(|2012al ). Asymmetries in bipolar outflows would have 
clear observational consequences as well, and the behav- 
ior of these outflows in different cooling regimes is an 
avenue of potential research. In any case future studies 
of the global structure of accretion flows should consider 
improvements such as fully 3D simulations, general rel- 
ativistic corrections, self-gravity, and more detailed heat 
and radiation transport mechanisms. One heat trans- 
port mechanism we have studied in some detail at low 
Eddington ratios is conducti on due to thermal electrons 
(jlnogamov fe Sunvaev Il2~010f h At lower Prandtl number 
convective overturning of hot gas can be suppressed, and 
the strength of the inflow and outflow can vary by a factor 
of a few. The distribution of outflow between polar and 
equatorial outflow also depends on the strength of the 
conduction, but the essential feature of mass accretion 
rates nearly two orders of magnitude below the Bondi 
rate remains unchanged. 

5. CONCLUSIONS 

Our high inflow rate solut ions are very similar to the 
standard and widely adopted IShakura fc Sunvaevl (|1973f ) 
results. But our low inflow results are to zeroth order 
the stationarv lPapaloizou fc Pringle I (|1984l ) solution. To 
next order in the small assumed viscosity they show pat- 
terns of circulation, with outflow almost balancing in- 
flow and the net accretion rate through a geometrically 
thick disk falling to only the order of a (the dimension- 
less viscous parameter) times the inflow rate. We label 
this behavior a RRIOS solution for "Radiating, Rotating 
Inflow-Outflow Solution" , as it is a significant generaliza- 
ti on of the ADIOS (Adiabatic Inflow-Outflow Solution 
of Blandfo rd fc Begelmanl (|1999[ )) applicable to the com- 
mon case wherein the angular momentum load of the 
inflowing gas strongly restricts the level of net accretion, 
at low inflow rates. In this low inflow state, viscous forces 
can drive an equatorial outflow that propagates out be- 
yond the Bondi radius. Further, entropy released by vis- 
cous stresses can accumulate near the centrifugal barrier 
and lead to episodes of convective overturning. Observa- 
tionally, the solutions that we find appear to be consis- 
tent with the fact that most nearby massive black holes 
show very low bolometric luminosity, considerably below 
that anticipated for spherically symmetric (non-rotating) 
Bondi-Hoyle flows, but also show occasional outbursts in- 
dicating a much higher rate of episodic outflow. 
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TABLE 1 

Resolution dependence of mass inflow rate at Rb, 

EQUATORIAL AND POLAR OUTFLOW RATES AT Rb, AND THE 
ACCRETION RATE FOR M B /M Edd = 10~ 4 . 



Resolution 


M ia /M B 


M polar /M B 


M cq /A/ B 


M acc /M B 


.81 


0.32 


0.034 


0.29 


0.0041 


1. 


0.27 


0.021 


0.24 


0.0042 


1.25 


0.27 


0.020 


0.25 


0.0044 


1.5 


0.24 


0.020 


0.22 


0.0045 


2 


0.21 


0.015 


0.18 


0.0040 



TABLE 2 

Resolution dependence of mass inflow rate at R b . 

EQUATORIAL AND POLAR OUTFLOW RATES AT Rb, AND THE 
ACCRETION RATE FOR M B /M Edd = 10 _1 . 



Resolution 


M in /M B 


M poUl /M B 


M cq /M B 


Ma.cc/ Mb 


.81 


0.51 


0.061 


0.004 


0.59 


1. 


0.45 


0.034 


0.014 


0.35 


1.25 


0.49 


0.033 


0.006 


0.40 


1.5 


0.55 


0.O40 


0.001 


0.48 


2 


0.59 


0.044 


0.001 


0.47 



APPENDIX 



In this Appendix wc provide further detail on the general properties of our solutions, including convergence tests. 
Figure H] illustrates the time dependence of characteristic quantities in our solutions for Eddington ratio Mb / ^Edd = 
10~ 4 and computed at our standard resolutions. We show the mass inflow rate at Rb, equatorial and polar outflow rates 
at Rb , and the accretion rate. The mass inflow rate, equatorial outflow rate, and accretion rate are statistically-steady 
over the illustrated interval. The polar outflow at large radii is punctuated by bursts of con vective overturning with 
high mass outflow r ates, but is also quasi-steady after time t ~ 55ie- The simulations of, e.g.. lStone et al. 1 (|1999[ ) and 
I Yuan et~ . (2012aj), found persistent convective eddies disrupting mass inflow near the centrifugal support radius, but 
we find that the thermal energy in the subsonic hot accretion flow with enforced symmetry through the equator can 
only support intermittent episodes of overturning extending to the Bondi radius. Table [T] shows the dependence of the 
four quasi-steady mass fluxes from Figure [9] on the spatial resolution of our simulations for the same Eddington ratio 
M B /M Edd = 1CT 4 . The first column gives the factor by which the number of grid points in each spatial dimension is 
multiplied. The accretion rate in standard resolution simulation is within 5% of the value at double the grid resolution, 
and the other mass fluxes are within ~ 30%. 

We perform the same exercise at Eddington ratio M B /M Edd = 10 _1 and show the quasi-steady average mass fluxes 
in Table El The mass inflow rate at Rb , polar outflow rate at Rb , and the accretion rate for the standard resolution 
simulation are within ~ 25% of the values in the simulation with double the grid resolution. Note that mass fluxes 
into and out of our grid do not have to exactly balance during a quasi-steady phase, as there can be a net buildup of 
gas near the centrifugal radius. This excess gas will eventually be ejected outwards during a convective overturning 
episode. 

We also explored the dependence of our solutions on our choice of no torque inner boundary condition located at 
R m = 91i? s . We have determined that mass fluxes at large radii change by a factor of 1.1, mass fluxes in the polar 
region at large radii increase by a factor of 1.2, and accretion rates decrease by a factor of 1.45, as we move the inner 
boundary to R- m = 91i? s /2 for solutions with Me/MEdd = 10~ 4 . Moving the inner boundary inwards by a factor 
of 4 to i?j n = 91R s /4 for the same Eddington ratio reduces the accretion rate by a factor of 2. There is ce rtainly a 
dependence of accretion flow eigenvalues on the location of the inner boundary (McKinncy fc Gammidl2002l ). but the 
differences do not fundamentally modify the solution types. This scal ing of accret i on rate wit h the location of th e 
inner and outer boundaries is also roughly consistent with that found by I Yuan et al. I (|2012aD and lStone et al. I (|1999l) . 
The transition to the cold disks is also pushed to a higher Eddington ratio as the inner boundary is moved inwards. 
A more self-consistent approach to determining the exact transition Eddington ratio would follow infalling gas down 
to the black hole horizon. 

Anpthcr variable that can affect the properties of the hydrodynamic accretion flow is the kinematic viscosity. We 
have explored other values of v for the hot disk solution at Eddington ratio Mb/M-eaa = 10~ 4 , and find that the mass 
inflow rate from large radii is proportional to y/v. The accretion rate as a fraction of the mass inflow rate from large 
radii appears to scale linearly with the kinematic viscosity. 

Finally, in Figure [TU] we show the scaling of density with radius for the hot disk solution at Ms/MEdd = 10~ 4 . 
Densities have been averaged over times t = 55ts — 90ts, and the solid curve represents a slice along the equator 
at 9 = it/2, whereas the dashed curve represents a slice at constant = 7r/4. The equatorial density scales roughly 
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as p oc r in the inner parts of the flow interior to the centrifugal support radius at R c = 0.02Rb- Beyond the 
centrifugal support radius, R c , the density begins falling more steeply as p oc r -12 , before levelling off in accordance 
with our outer boundary condition. In comparison, the density scales as p oc i?~ 3 / 2 in ADAF models (jNaravan fc Yil 
I1994T ) and as p oc R^ 1 / 2 in standard CDAF models (jQuataert fc Gruzinovll200Cll : iNaravan et al.ll2000f ). though we use 
a different viscosity prescription. 
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Fig. 9. — Time dependence of mass inflow rate at Rb, equatorial and polar outflow rates at Rb, and the accretion rate for the solution 
at Mfl/^Edd = 10 -4 . The simulation has converged to a unique steady-state solution after t ~ 55te- 
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